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ABSTRACT 

Theoretical investigations have suggested the presence of Intermediate Mass 
Black Holes (IMBHs, with masses in the 100 — 10000M Q range) in the cores of 
some Globular Clusters (GCs). In this paper we present the first application 
of a new technique to determine the presence or absence of a central IMBH in 
globular clusters that have reached energy equipartition via two-body relaxation. 
The method is based on the measurement of the radial profile for the average 
mass of stars in the system, using the fact that a quenching of mass segregation 
is expected when an IMBH is present. Here we measure the radial profile of mass 
segregation using main-sequence stars for the globular cluster NGC 2298 from 
resolved source photometry based on HST-ACS data. NGC 2298 is one of the 
smallest galactic globular clusters, thus not only it is dynamically relaxed but also 
a single ACS field of view extends to about twice its half-light radius, providing 
optimal radial coverage. The observations are compared to expectations from 
direct N-body simulations of the dynamics of star clusters with and without an 
IMBH. The mass segregation profile for NGC 2298 is quantitatively matched 
to that inferred from simulations without a central massive object over all the 
radial range probed by the observations, that is from the center to about two 
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half-mass radii. Profiles from simulations containing an IMBH more massive 
than « 300 — 5OOM (depending on the assumed total mass of NGC 2298) are 
instead inconsistent with the data at about 3a confidence, irrespective of the 
IMF and binary fraction chosen for these runs. Our finding is consistent with the 
currently favored formation scenarios for IMBHs in GCs, which are not likely to 
apply to NGC2298 due to its modest total mass. While providing a null result 
in the quest of detecting a central black hole in globular clusters, the data-model 
comparison carried out here demonstrates the feasibility of the method which 
can also be applied to other globular clusters with resolved photometry in their 
cores. 

Subject headings: stellar dynamics - - globular clusters: general - - methods: 
n-body simulations 



Introduction 



Intermediate Mass Black Holes (IMBHs — with masses of order 1 2 M M to 10 4 M^) have 



been suggested to form in the cores of young star clusters (e .g., see iPortegies Zwart et al. 
2004 ) and in dense globular clusters (IMiller &: Hamilton! |2002| ). Possible observable features 
related to these objects such as shallow cusps in the density and velo city dispersion profi le of 
the cluster stars were proposed early on (e.g . see the seminal paper of lFrank fc Reedll976l and 
van der Marell 12004] ; IMiller fc Colbertll2004l for a review). Cuspy central Surface Brightness 
Profiles (SBPs) have now actually been observed in HST data on a sizable fraction of GCs 
( jNoyola fc Gebhardtll2006l ). but it is unclear if they are directly related to IMBHs. Detections 
based on the analysis of line-of-sight velocity data have been made but the ev idence does not 
appear to be conclusive (see the introductory discussion in I Gill et al.l 120081 ). For example, 
a paradigmatic case is that o 
extragalactic cluster Gl (see 



the much debated detection of a 2 ■ 10 4 M t r ) black hole in the 



Gebhardt et al 



20051 : iTrudolyubov &: Priedhorsky 



2004 



2002 



Ulvestad et al. 



2003). 



Baumgardt et al.ll2003l ; iGebhardt et al. 



Ideally a direct unambiguous detection of an IMBH is possible in GCs by measuring 
orbits of stars bound to the BH, which would also allow a precise measurement of the 
central mass. Present-day HST imaging capability has the accuracy required for this kind of 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. This paper is associated with program #11284. 
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observation but a significant investment of time spread over multiple epochs is required. It 
is therefore necessary to find preliminary criteria to narrow down the list of candidate GCs 
for focused follow-up observations. 

In this paper we apply to NGC 2298 a ne w method we recently proposed for assessing 
the presence of an IMBH (see iGill et al.ll2008l ). The idea is to quantify the amount of mass 
segregation present in a collisionally well-relaxed stellar system, that is with a half-light 
two-body relaxation time below one billion years. We have shown through direct N-body 
simulations that the presence of an IMBH heavily affects mass segregation of stars in a GC. 
Systems hosting an IMBH develop a low degree of mass segregation, as opposed to IMBH-free 
GCs in which more massive stars move preferentially towards the center of the cluster over a 
relaxation timescale. The presence of the IMBH tends to equalize the velocity dispersions of 
all stellar mass components in the system, thus reducing radial mass segregation. Note that 
the method we propose is not applicable to the most massive galactic globular clusters such 
as Gl and Q Centauri, because their two-body relaxation time is too long and the amount of 
observed mass segregation does not necessarily reflect its long term equilibrium value. Our 
investigation thus aims to cover a different region of galactic GC par ameter space compared 
to t he one explored by c urrent claims of GCs IMBH detection such as lGebhardt et al.l (120021 ) 
and iNovola et all £53). 



The radial profile for mass segregation in main-sequence (hereafter MS) stars can be 
readily measured in GCs via observations with sufficie: 
ual sources in the crowded cores of these systems (e.g. see 



De Marchi & Pulone! 


20071). In this 


as analyzed by 


De Marchi & Pulone 



( 120071 ). These data are ideal to measure the radial variations of mass segregation because 
they range from the center out to more than twice the half-light radius. The observed radial 
variation in t he mean MS mas s is then compared to expectations based on the numerical 



simulations of IGill et al.l (120081 ). which allow us to constrain the mass of a central BH in the 
system. 

This paper is organized as follows. In Section [2] we discuss the HST observations for 
NGC 2298 and our data analysis, which is then compared to our N-body simulations in 
Section [3J We summarize and conclude in Section HI 



2. Observations and Data Analysis 



We study the mass segregation of MS stars of NGC 2298 using deep HST/ ACS ob- 
servations in the F606W and F814W bands. The field covers an area of 3.4 • 3.4 arcmin 2 
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around the c luster center, extending to more than t wice the cluster's hal f- light radius (see 



Harris! Il996l ). The data reduction is described in |Pe Marchi &: Pulond (120071 ). who de- 
rive the color-magnitude diagram that we use in our analysis. The data have 10<r de- 
tection limits of m F606W ~ 26.5 and uifsuw ~ 25 with a completeness at the detection 
limi t above 50%. Background sources contamination is essentially negligible for this clus- 
ter (IDe Marchi fc Pulond 120071 ). We assume a distance m odulus of 15. 15 mag (distance 
of about 12.6 kpc), and color excess E(B — V) = 0.14 from lHarrisI (119961 ). From the color 



magnitude diagram we in 
luminosity relat ion fro m 



t fer the mass of each i ndividual main sequence source using the mass- 
Bar affe et al.l (119971 ) and assuming a metallicity [Fe/H] = —1.85 
for the cluster (jHarria 119961 ) . Our MS stars catalog consists of objects with masses in the 
range [0.2 : O.8]M , where the lower limit is set by completeness cutoff and the upper limit 
by the turn-off mass. 



2.1. Cluster Properties 



The Surface Brightness Profile for NGC 2298 is part of the iTrager et al.l (119951) com- 
pilation of photometric data on galactic GCs. McLaughlin fc van der Marel 12005) fitted a 



De Marchi fe Pulone 



single -mass King model to it, obtaining a total mass of 3.O9-1O 4 M while^_ 

( 2007 ) derive a total mass of 5 • 1O 4 M based on a multi-mass dynamic model. The lHarrisI 

( 1996 ) catalog reports a projected half-light radius of 46.8 arc sec. while iMcLaughlin fc van der Marel 
( 2005 ) derive ru = 45.4 arcsec, corresponding to 2.35 parsecs. De Marchi fc Pulone! ( 20071 ) 
find instead a larger half-mass radius — that is rh m = 72 arcsec — from their multi-mass 
model. 

From the structural par ameters of NGC2 298 we can derive its half-light relaxation time, 
defined in physical units as (lDjorgovskilll993l ): 



' rti 



8.9- 10 V / 1M G 



x 



log(0.4iV) \{m # ) 



Mi 



tot 



1M 



© 



0.5 



r h i 
lpc 



1.5 



where (m*) = M tot /N is the average mass of a star (including dark remnants), iV is the 
number of stars and M tot is the total mass of the system. 

Assuming an averag e stellar mass of 0.5M^, the half-lig ht relaxation time for NGC 
2298 is 10 8 41 yr usin g the McLaughlin fc van der Marell (120051 ) structural parameters. If we 
consider instead the De Marchi fc Pulone ( 20071 ) modeling R Eq. CE] yields t T h = 10 8 76 yr. 



2 Note that the relaxation time of 3.1 Gyr quoted in be Marchi fc Pulond (|2007h for NGC 2298 is the 
result of a typo in the output of the code they used to contruct the dynamic model. 
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In both cases NGC2298 appears to be a dynamically old cluster even af ter accounting for 



a sign ificant mass loss in the cluster during its evolution, as suggested by IBaumgardt et al. 
( 2008 ): if the cluster was originally four times as massive and twice as large there would 



still have been enough time to relax, as the initial relaxation time would have been t r / l (t = 
0) ~ 3Gyr. 

Further evidence for an old dynamic age comes from the shape of the mass function 
of NGC 2298, which shows a depletion of low mass stars. This makes the cluster a good 
candidate for our search because it would be expected to have undergone core collapse, unless 
there is a central source of energy capable of halting the core contraction, such as an IMBH 

l2007bO l. While 



Trenti et al 



or a significant population of primordial binaries (e.g. see 
the current formation scenarios for IMBHs in GCs assume a higher clust er mass than the 
curre nt NGC2298 mass, its probable initial mass — up to 1.4 x 1O 5 M (IBaumgardt et al. 
20081 ) — was well in the range required for IMBH formation. 



In order to consistently compare the data with our sample of numerical simulations, 
we identify the projected half-mass radius in main sequence stars within the ACS field of 
view — which we call Th m . Given that the mass-to-light ratio of globular clusters has radial 
gradients, this quantity is different from the standard half-light radius used in the literature. 
To measure r^ m we bin the star counts in cluster-centric radius and sum over star mass, 
applying the completeness correction appropriate for the given cluster-centric radius and star 
magnitude. We obtain the surface density profile of the cluster from mai n sequence stars and 
then a pply the non-parametric spline-smoothing technique described in iPasquato fc Bertin 
(120081 ) to the profile obtaining our best estimate of the total main sequence mass of the 
cluster in the field of view and of the respective half-mass radius. We do not extrapolate the 
light profile outside the ACS field which might lead to the discrepancies in the determination 
of the structure of the cluster discussed above, but rather include the effect of a finite field 
of view into the simulation analys i s (see Section [3]). We derive r hm = 49.0 arcsec, slightly 
larger than ru reported in lHarrisI (119961 ) . This is not surprising, as the light profile tends 
to be dominated by red giant stars, more massive on ay erage than MS stars and therefore 
more centrally segregated (for example, see iHurleyi 120071 ) . 



Trenti et al.l (|2007bf ) find that the half-mass radius stays approximately constant while a star cluster is 



being tidally disrupted 
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2.2. Mass Segregation 



We quantify mass segregation through the radial variation of the mean mass {w) MS 
of MS stars. At the center of the system, where high mass stars segregate through energy 
equipartition, the mean m ass of stars i n such range is expected to be higher than in the 
outer region of the system (ISpitzerlll9871 ). To construct (Tn) MS (r) we apply the completeness 
correction as follows: 



(m) M s(r) 



(2) 



J2ii/f(m,ri) 

where Mj is the star's mass and f(m,r) is the completeness for a star of magnitude m 
at cluster-centric 2D radius r. To apply the completeness correction, we assumed a one 
to one map from mass to luminosity (Mj to for main-sequence stars (specifically that 
of iBaraffe et al.lll997l ). The color information is of course used to select main sequence 
stars. The sum over the index i is carried out over all stars in an annulus around the 
radius r. Since the complet eness function was defined in both radial and magnitude bins by 
De Marchi &: Pulond (120071 ). we define a continuous f(m, r) by bilinear interpolation. As a 
compromise between obtaining high spatial resolution and minimizing Poisson fluctuations 



we calculate (m) 



using 20 concentric annuli comprising each 5% of the stars in number. 



Errors on each radial annulus are obtained using a bootstrap technique (lEfroru Il979l ) . 
We proceed as follows: Let {Mj} i=1A r be the catalog of observed main sequence star masses 
within a radial annulus. We use this catalog to generate 100 synthetic catalogs. To construct 
a synthetic catalog, we first extract N random numbers {q^j^^N from a uniform random 
number generator in the range [1 : N]. The synthetic catalog of masses is then defined as 
{M q .}j = i^, that is we extracted with replacement and uniform probability from the original 
data-set. Eq. [2] is then applied to each synthetic catalog of a radial annulus and the 1 a 
error on the observed {m) MS (r) is defined as the standard deviation of the sample of the 100 
synthetic catalogs for that radial position. 

The observed mass segregation profile, normalized to the average mass measured around 
Thm (using the average MS mass measured between 0.8 and 1.2 r^ m ) is indicated as Ara(r) 
and shown in Fig. [IJ From the plot it is clear that the cluster has a marked radial variation 
in the MS mass, with a difference of about 0.14M Q from the center to the outermost region 
at r > 2r/ MTi . This is about 30% of the average mean MS mass of the system O.529±O.OO2M . 
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3. Comparison with Numerical Simulations 



The numerical simulations used in this paper have been carried out wit h a state-of-the- 
art direct N-body code for sta r cluster dynamics, NBODY6 (jAarsethll2003l ). NBODY6 has 
been modified as discussed in iTrenti et al.l (j2007al ) to improve accuracy in the presence of 
an IMBH, and uses regularization of close gravitational encounters without any softening, 
guaranteeing high accuracy on the integration of the star trajectories even in the proximity 
of the black hole. 

The det ails of the N-bod y simulations — along with a mass segregation analysis — are 
discussed in iGill et al.l ( 120081 ). here we summarize their properties (see also Tab. HJ). We 
use from 16384 to 327 68 particles. T he ir individual masses are drawn from an initial mass 



function (IMF, either ISalpeterl Il955l or iMiller fc Scald Il979l ) and then evolved thro u gh an 



instantaneous step of stellar evolution to ~ 12 Gyr of age using the iHurley et al.l (120001 ) 
tracks. Our stellar mass black holes are in the [5 : 10]M o range. The stellar evolution step 
evolves the cluster IMF to match the turn-off mass of ~ 0.8M Q observed in NGC 2298. 
Simulations also have up to 10% of primordial binaries and include the effects of a galactic 
tidal field. 

Simulations are carried out until complete tidal disruption of the systems, which happens 
after about 20 initial half-light relaxation times t r h (defined in Eq.[T]), a timescale l onger than 
the age of the Universe for a typical globular cluster. As discussed in I Gill et al.l (120081 ). the 
simulations settle down in a quasi-equilibrium configuration with respect to mass segregation 
after a few relaxation times. At this point, there are three regimes for the asymptotic value 
of mass segregation measured as the difference in the main sequence mass between the 
center and the half mass radius^ (A(m)): a low value (A(m) < 0.07) corresponding to 
an high probability of harboring an IMBH, a high value (A(m) > 0.1) associated to a 
low probability of harboring an IMBH and an intermediate regime, where models with and 
without an IMBH both can be found. In addition, the mass segregation profile measured 
through Am(r) and/or A(m) is a differential measure normalized to the half-mass radius, 
so the averag e value of the M S star mass over the whole cluster is not important and runs 
with either a lSalpeterlll955l or lMiller fc Scalolll979l IMF have similar Am(r) profiles. 



Furthermore, our N=32k runs starting with a lMiller Sz Scald (119791 ) mass function have 



4 Note that the quantity A(m) defined in IGill et al.l (|2008f ) has a slightly different normalization with 
respect to Am(r = 0) used in this work. We define Am(r = 0) as th e difference in me an MS mass between 
the center and mean MS mass in the radial annulus r/r/j m £ [0.8, 1.2]. I Gill et al.1 (|2008l ) use instead a smaller 
annulus around the half mass radius, defined as that containing 5% of mass of the cluster and centered on 
Thm- In Table 1 we report both indicators for comparison. 
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a similar shape of the mass function in mai n sequence once they ha ve lost about 75% of their 
initial mass (as estimated for NGC 2298 by iBaumgardt et al.ll2008l ). This is shown in the left 
panel of Fig. [2]for our simulation 32km at t — lQt r h- dynamic evolution of the cluster and the 
subsequent preferential ejection of light stars has naturally evolved the initial Miller & Scalo 
distribution to match that observed in NGC 2298, especially for the simulation with no BH. 
To ensure a self-consistent comparison of the observations with our simulation, Fig. [2] has 
been obtained by taking into account the completeness limit of the data when "observing" 
the simulation: we have projected the simulation snapshots along random lines of sight, then 
rejected with probability 1 — f(m,r) particles with magnitude m at a 2D radius r. 

To construct the mass segregation profile of our simulations, we use all the snapshots 
from lGill et al.l (120081 ) with N > 16384 particles. We restrict the analysis to the time interval 
between 7 and 9t r h, a period in which the simulated clusters have reached their asymptotic 
amount of mass segregation but still have more than half of the initial number of particles 
(except for the plot shown in Fig. [2J obtained using snapshots between 15.5 and 16.5 t rh ). 
This gives us 26 snapshots per simulation for runs with N = 16384, and 32 snapshots for 
runs with N = 32768 for a total of 324 snapshots available for comparison with observations. 
As in this paper we aim at analyzing the simulations as similarly as possible to the NGC 
2298 data, we recast the analysis in terms of the full radial profile for mass segregation 
rather than limiting to A(m) as in lGill et al.l (120081 ). Only line-of-sight projected quantities 
are used and only main sequence stars are considered in the analysis. For those runs with 
primordial binaries, binary stars are projected onto the main sequence for single stars and 
treated as a sin gle star with mass equal to that of the heavier main sequence member (see 
Gill et al.ll2008l ). The analysis is also restricted to MS stars more massive than 0.2M Q to 
match the completeness limit of the observations. Under these assumptions we calculate 
the center of visible projected mass for each simulation snapshot (this turns out not to 
differ significantly from the center of mass of all gravitationally bound particles) and a first 
guess of r hm . We then repeat our analysis including only particles lying within 2r hm from 
the center and recompute both the center of visible mass and the final r^ m value. This 
procedure closely resembles the observational limit on the field of view of our data, which 
extends to about 2r^ m . To bring the number of particles "observed" in a simulation snapshot 
close to the actual number of stars in NGC 2298 we sum three independent projections for 
each snapshot. In this way we almost reach a 1:1 ratio of particles to stars for our larger 
N = 32768 runs. 



Fig. [T] shows the results from the lGill et al.l (120081 ) sample of snapshots superimposed to 
the observed mass segregation profile derived in Section [2J The blue shaded area corresponds 
to the region defined by the 2a contours for points derived from all the snapshots of all runs 
with an IMBH (thus encompassing simulations with different IMF and binary fractions as 
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reported in Table 1), while the green shaded area is the 2a region for runs with no central 
BH (including all the no BH runs in Table 1). 

The data-points for NGC 2298 are fully encompassed by profiles derived by simulations 
run without a black hole and the central profile (innermost two measurements) is outside 
the region associated to simulations with a central BH. The la error bars associated to the 
observed points are smaller than the scatter from simulation snapshots, thus the main factor 
in setting the confidence level at which we can exclude a central black-hole is set by the 
variance in the simulations. The upper right inset of Fig. [1] shows the upper envelope of all 
162 mass segregation profiles obtained from snapshots of simulations containing an IMBH. 
The two central observed data points lie above such envelope at a combined confidence 
level greater than 2a confidence level, providing us with a quite stringent test against the 
possibility that the observed high amount of mass segregation in NGC 2298 is the result 
of a random fluctuation of a system with a central IMBH. If this were the case, two > 2a 
fluctuations both in the simulations and in the observations would be required, thus we can 
exclude this scenario at about 3a level. 

To better quantify the influence of different initial conditions on the expected variance of 
the mass segregation profile, we used a second, different approach. Rather than calculating 
the variance of all sample of simulations with or without a central IMBH, we compute the 
variance in Am(r) separately for each run. In Fig. [3] we then define the blue and green 
shaded areas as the envelope of the 2a regions of each run. In this way simulations with 
different IMF and binary fractions are treated separately and the systematic differences 
between them are not treated as random error. The two shaded areas are therefore larger 
by construction than in Fig. [TJ Still a clear separation is present between runs with and 
without an IMBH. showing that systematics due to different IMF and binary fractions are 
not a concern for our results. 



Fig. [3] is also based on four additional runs with respect to the lGill et al.l (120081 ) runs that 
were used to construct Fig. [IJ These additional runs, listed a t the end o f Table 1, have been 
generated using the same initial conditions of run 16kmbh in iGill et al.l (120081 ) . except for a 
different seed of the random number generator. These runs have been added to further study 
mass segre gation in p r esence of an IMBH when there is a Miller & Scalo IMF. Note that run 
16kmbh in IGill et al.l (120081 ) is the one which shows the highest degree of mass segregation 
among th o se wi th an IMBH, even though both its N = 8k and iV = 32k counterparts in 
Gill et al.l (120081 ) have a lower degree mass segregation. Our additional iV = 16k runs show 
on average a slightly lower degree of mass segregation as report ed in Tabl e 1 (eig ht column), 
even though one of the four runs is very close to the original IGill et al.l (120081 ) results. In 
preparing Fig. [3] we treated all these five 16kmbh runs independently, thus the 2a envelope 
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shown in the figure includes a > 99% confidence level area with respect to this specific 
initial condition. Note also that the addition of 106 snapshots from four new runs have 
naturally lead to a slight increase of the maximum amount of mass segregation measured in 
the simulations (plotted in the upper right inset of Fig. [3D compared to the inset of Fig. [TJ 
But these additional snapshots also increase the confidence level for BH rejection: the inset 
in Fig. [3] is based on a total of 266 snapshots and still none of those snapshots reaches the 
data whithin the observed error bar (in particular with respect to the second innermost data 
point). 

If we restrict our analysis to snapshots that have a global mass function similar to 
that observed in NGC 2298, our results are stronger. The right panel of Fig. [2] shows a 
comparison of the observed AM(r) with the expectations from our N=32k Miller&Scalo 
simulations evolved around t = 16t r ^, when iV ~ 8000 parti cles remains in the syst em (a 



mass loss consistent with the one estimated for NGC 2298 by lBaumgardt et al.l 120081 ) . The 
global mass function of the simulations at that point is an excellent match to the observed 
one, especially for the run without central BH. Evaporation of a significant fraction of stars 
has marginally lowered the amount of central mass segregation compared to that measured 
at t m 8t r h for the same runs. The contour areas for the simulated AM(r) in this figure are 
contained within those of Fig. [TJ but are significantly more compact in the central region. 
Therefore from Fig. [2] the presence of a central BH can be rejected at a significance level much 
higher than 3cr. In addition, we get a fully consistent representation of the mass segregation 
profile from the simulations without a central BH. The quantitative match of the observed 
mass segregation profile to the numerical expectations also provide a posteriori evidence that 
NGC 2298 has reached its equilibrium value for mass segregation and thus it is collisionally 
relaxed (see discussion in section |2~TT]) . 

All our simulations that include an IMBH have a ratio of about 1% between the total 
cluster mass and the black hole mass. Assuming a total mass for NGC 2298 of 3.09 x 10 4 M Q , 
this implies that we can exclude the presence of a central BH of mass > 3OOM . The latter 
limit would increase to 5OOM if the total mass of NGC2298 is instead ?s5x 1O 4 M . 



3.1. Systematic Uncertainties 

In order to compare simulations to observations there is need to fix a radial scaling. 
We choose to normalize the mass segregation profile Am(r) to its value at r^ m , as discussed 
in Section [2] and shown in Fig. [TJ A correct measurement of r^ m is therefore vital for a 
consistent comparison between simulations and observations. For the simulations, we have a 
complete control over the observables in the system, but we have to evaluate possible biases 
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in the measure of r^ m from the photometry of NGC 2298. Mis-estimating the value of mean 
MS mass at r hm due either to local fluctuations of (m) MS (r hm ) or to a wrong estimate of r hm 
can then in principle lead to an un-physical shift of the whole profile upwards or downwards 
along the vertical axis. 

Local fluctuations in the average main sequence mass are reduced by our choice to 
normalize the Am(r) profile using the mean MS mass between 0.8r/j m and 1.2r/i m , instead 
of the punctual value at exactly r hm (see Section [2]). We note in passing that also choosing 
a smaller interval ([0.9, l.ljr^) yields comparable results, although with more noise. 

To assess the impact of a mis-determination of r^ m , we repeat our analysis using values 
for rh m shifted by ±4 arcsec from our best estimate r^ m = 49. A change in rh m of this 
amount does not critically affect our conclusion, as it shifts Am by less than the la error 
associated to the measure. The amount of observed mass segregation is decreased if one 
adopts a smaller value for r^ m . Only a significant change in rhm ijhm % 40 arcsec) would 
move the measurement within the amount of segregation typically associated to the presence 
on an IMBH. 

Another source of syst ematic uncertainty that nee ds to be addressed is the determination 



of the center of the cluster. iNoyola &: Gebhardtl (120061 ) point out that center determination of 



GCs is a difficult problem, with literature/ground based coordinates of GC centers sometimes 
being inaccurate. Miscentering of the cluster in the analysis of observational data could in 
principle lead to an artificially shallower mass segregation profile in the central part of the 
cluster. If this is the case, then the confidence level of our null result would just increase. 
However, miscentering is not likely to be a significant issue in our case. Not only did we 
use high resolution data from HST but we determined the center based on the mass of main 
sequence stars, not on the total light from the cluster. Therefore the Poisson fluctuations 
in the small number of red giant stars that might dominate the light profile of a globular 
cluster do not affect significantly our analysis. To quantify the random error on our adopted 
center position we carried out a Monte Carlo bootstrap resampling test. We generated 100 
synthetic samples of stars by extracting with replacement main sequence stars from the 
observed catalog. We then re-computed the center of each synthetic sample and used the 
100 center coordinates thus obtained to calculate the la fluctuation of the center position. 
The value we recover from our Monte Carlo test is below 0.4 arcsec. 



-12- 



3.2. Applicability of mass segregation indicator to constrain IMBH presence 



In iGill et al.l (120081 ) we suggest two conservative criteria to ensure that a stellar system 
has reached its long-term mass segregation profile: (1) have a tidal to half-light radius 
r t/ r h ^ 10 and (2) have an half-light relaxation time below 1 Gyr. NGC 2298 fulfills the 
second condition, even assuming a worst case estimation of its structural para meters (see 
Sectio n [2TTT) . but marginally fails the first one, as rt/rh = 8.3. As discussed in I Gill et al. 



( 120081 ). these requirements are sufficient, but not necessary for the mass segregation to 
have reached its equilibrium value. Based on a detailed dynamic modeling of the cluster 
in Section 12.11 we demonstrate that the cluster appears to be well relaxed and its stars to 
have reached energy equipartition, despite its significant mass loss due to tidal evaporation. 
Furthermore, if there is some primor dial mass segregation, like it has been suggested for 
NGC 2298 by lBaumgardt et al.l (120081 ). then the time needed to develop the asymptotic mass 
segregation profile becomes shorter. This provides a further validation for the application of 
our analysis to this system. 



4. Discussions and Conclusions 

We have analyzed ACS-HST data covering the globular cluster NGC 2298 to quantify 
the radial variation of mass segregation, measured from main sequence stars in the range 
[0.2 : O.8]M . For dynamically relaxed systems such as NGC 2298 (that is with a half-mass 
two-body relaxation time well below the Hubble time), the degree of mass se gregation presen t 



P 

1. : 



can shed light on the presence of a central intermediate mass black hole (see lGill et al.ll2008l ). 
The observed mass segregation profile (see Fig. [1]) has been compared to the expectation from 
a sample of direct N-body simulations with and without a central IMBH. The simulations 
have been analyzed as closely as possible to the observational data. We find that NGC 2298 
exhibits a fair amount of mass segregation, with the average main sequence stellar mass be- 
ing larger at the center compared to that at the half-mass radius by Am(r = 0) ~ O.O95M . 
This observed radial variation of the main sequence mass is fully encompassed by curves 
derived from simulations without a black hole. None of the simulation snapshots with a 
central IMBH presents such a large degree of mass segregation and the two innermost data- 
points lie each at 2a above the maximum mass segregation measured in the simulations, 
for a nominal combined confidence level of about 3a. At this confidence level we can thus 
consider unlikely the presence of an IMBH of mass Mbh ^ 300 — 500M Q (depending on 
the total mass of the system). Certainly NGC 2298 does not appear a promising candidate 
for follow-up observations to search for an IMBH such as with a proper motion study. This 
could partly have been expected based on the small size of NGC 2298, because the proposed 



13 



scenarios for IMBH formati on in GCs require a higher cluster mass (jMiller fe Hamirtonll2002 
Portegies Zwart et al.ll2004l ) than its present-day mass. On the other hand, dynamical inter- 
action with t he galactic environment is likely to have stripped the majority of the cluster's 
initial mass (IBaumgardt et al.ll2008l ). which could have been well above the threshold for 
runaway collapse of massive stars. 

A critical assumption in our analysis is a proper evaluation of the current dynamic state 
of NGC2298 and in particular of its half-light relaxation time, which is complicated by the 
presence of potentially inconsistent structural parameters in the literature (see our discus- 
sion in Section [2. II) . Even assuming a worst case scenario with respect to the structural 
parameters published, the current half-light relaxation time we derive is below ~ 0.6Gyr, 
guaranteeing that the cluster had enough time to reach energy equipartition and its equi- 
librium value of mass segregation, even if it has lost most of its initial mass and thus its 
initial half-ligh t relaxation time could have been as high as ~ 3 Gyr. With respect to this 
issue, note that iGill et al.l (120081 ) recommend to apply the mass segregation analysis to well 
relaxed clusters, defined as having a current t r h < lGyr and r t /rhi > 10 in order to avoid 
false positives, that is clusters that have a low amount of mass segregation not because there 
is an IMBH but because the cluster has not yet fully developed mass segregation. In the case 
of NGC2289 we are in the opposite regime as we find more mass segregation than expected 
in presence of an IMBH. Furthermore, we have a quantitative match of the observed radial 
profile of mass segregation with the expectations from our numerical models. This is a fur- 
ther evidence that NGC2298 is collisionally relaxed and our analysis of the full radial pr ofile 
of mass segregation goes beyond the single point measure discussed in IGill et al.l ( 120081 ). 



The absence of a central IMBH in NGC2298 has an interesting consequence for the 
dynamics of the system. The cluster has a rather large c ore and does not appea r to have un- 
dergone core collapse despite its advanced dynamic age (IDe Marchi et al.ll2007l ). This means 
that there must be another source in the core capable of generating kinetic energy through 
gravitational enc ounters, such as a population of primor dial binaries with a number density 



> 0.05 (e.g. see iHeggie et al.l 120061 ; iTrenti et al.ll2007bl ). The color-magnitude diagram for 



the cluster has indeed a widened main sequence which allows us to get a lower limit of 0.04 
for the binary fraction in the core and a likely fraction two to three times higher. 

We have addressed possible sources of systematic errors in constructing the observed 
mass segregation profile, namely mis-estimation of the half-mass radius r^ m and cluster 
miscentering (see Section 13.11) . Neither of these is expected to have a significant impact 
on our measure. Quantifying the robustness of the measured run-to-run varia tions from 
the sa mple of our simulations is instead more challenging. While the original IGill et al. 



( 120081 ) simulations have 162 snapshots with an IMBH and it is thus hard to solidly assess 
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fluctuations above the 2a level, Fig. [3] shows that the inclusion of 104 more snapshots from 
runs with an IMBH does not alter our main results, thereby showing that both systematic 
and statistic errors are under control. 

Despite the null result of this search, we have demonstrated the practical application 
of mass segregation as a fingerprint for the presence or absence of an IMBH. The method 
can be applied to other galactic globular clusters provided that data of HST quality for the 
central regions of the systems are complemented by the acquisition of a field located around 
the half-light radius of the system. 

We thank Enrico Vesperini for useful discussions and suggestions and the referee for 
a thorough and constructive report. Support for proposal HST-AR-11284 was provided by 
NASA through a grant from STScI, which is operated by AURA, Inc., under NASA contract 
NAS 5-26555. This work was supported in part through NASA ATFP grant NNX08AH29G 
and by the National Science Foundation under Grant No. PHY05-51164. 
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Fig. 1. — Observed radial mass segregation profile (Ara(r) measured in M ) for NGC 2298 
(red points with la error bars), compared to expectations from numerical simulations. In the 
main panel the blue (IMBH) and green (NO BH) shaded areas represent the 2a confidence 
level area for th e profiles from the ensemble of snapshots with iV > 16384 published in 



Gill et al.l (120081 ) (see Table 1 for a description of these runs). Also shown as long dashed 



lines are the inner la regions. The small inset shows the inner observed data points compared 
against the upper envelope of all the profiles associated to snapshots with a central IMBH, 
that is against the maximum mass segregation measured in the simulations with an IMBH. 
Mass segregation in NGC 2298 appears typical for a system without a central BH. 
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Fig. 2. — Left panel: Observed main-sequence mass function for NGC 2298 (red line) com- 
pared to that of our 32km and 32kmbh (NO BH: green line, IMBH: blue line) simulations 
at time t = lQt r h wh e n abo ut 75 % of the initial mass is lost, as estimated for NGC 2298 
by iBaumgardt et al. (120081 ). The simulations started with N = 32768 particles from a 



Miller fc Scale! (Il979l ) IMF and have been projected in 2D and "observed" with a field of 
view extending to 2r hm assuming the completeness of the NGC 2298 data (see section [3]). 
The IMBH simulation has a mass function less depleted of light stars because of its reduced 
mass segregation. Right panel: Observed radial mass segregation profile (Am(r)) like in 
Fig. [1] but with theoretical expectations based only on the N=32k snapshots taken in the 
time interval t e [15.5 : 16.5] t r h- 
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Fig. 3. — Observed radial mass segregation profile for NGC 2298, as in Fig. [U The blue 
(IMBH) and green (NO BH) shaded areas are constructed here using a different approach 
compared to Fig. [p (i) we include four additional runs with a central IMBH (whose initial 
conditions lead in lGill et al. 2008 to the highest amount of observed mass segregation; see 
detail in Section. [3]); (ii) for each simulation run we first define its 2a confidence level area 
for mass segregation and then construct the global 2a areas for the ensembles with and 
without a central BH taking the envelope of the individual runs confidence regions. This 
procedure better highlights run to run variations in mass segregation associated to different 
initial conditions and by construction results in larger uncertainty regions. It is therefore 
an extremely conservative approach and represents a stronger test than Fig. [T] to reject 
the presence of a central BH. The upper small inset is defined as in Fig. but here 104 
additional snapshots from the 4 additional runs are considered: this resulted in an increase 
of the maximum amount of mass segregation seen in simulations, still none of the sna/nshots 
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Table 1. Summary of the N-body simulations. 



TD 


N 


TMF 


1V1 IMBH 1 IVI tot 


M IMBH /M Q 


Jb 


(Am) 


I Am) 


&P09 


Am min 
L±rri P09 




16ks 


16384 


Sal 


N/A 




N/A 





0.11 


0.101 


0.008 


0.083 


0.113 


16ks.l 


16384 


Sal 


N/A 




N/A 





0.14 


0.13 


0.01 


0.112 


0.158 


16km 


16384 


M&S 


N/A 




N/A 





0.14 


0.137 


0.009 


0.116 


0.151 


16kbs 


16384 


Sal 


N/A 




N/A 


0.1 


0.09 


0.074 


0.008 


0.056 


0.090 


16kbm 


16384 


M&S 


N/A 




N/A 


0.1 


0.10 


0.102 


0.009 


0.085 


0.124 


32km 


32768 


M&S 


N/A 




N/A 





0.14 


0.142 


0.007 


0.128 


0.154 


16ksbh 


16385 


Sal 


0.015 




103.1 





0.05 


0.048 


0.006 


0.037 


0.060 


16ksbh.l 


16385 


Sal 


0.015 




60.9 





0.06 


0.06 


0.01 


0.041 


0.078 


16kmbh 


16385 


M&S 


0.015 




128.2 





0.08 


0.071 


0.009 


0.048 


0.091 


16kbsbh 


16385 


Sal 


0.01 




113.4 


0.1 


0.04 


0.042 


0.005 


0.033 


0.052 


16kbmbh 


16385 


M&S 


0.01 




141.0 


0.1 


0.05 


0.050 


0.008 


0.038 


0.072 


32kmbh 


32769 


M&S 


0.01 




240.0 





0.07 


0.069 


0.005 


0.058 


0.083 


16kmbh 7 


16385 


M&S 


0.015 




128.2 





N/A 


0.069 


0.008 


0.055 


0.083 


16kmbh 7/ 


16385 


M&S 


0.015 




128.2 





N/A 


0.063 


0.010 


0.040 


0.079 


16kmbh m 


16385 


M&S 


0.015 




128.2 





N/A 


0.055 


0.009 


0.042 


0.074 


16kmbh /y 


16385 


M&S 


0.015 




128.2 





N/A 


0.074 


0.009 


0.060 


0.093 



Note. - Properties of the N-body simulations used in this paper. The first column is the simulation ID 
(subscript .1 means that the IMF was dow n to O.IM^), the secon d column reports the number of particles in 
the run, the third the IMF (ISalpeterl Il955l or iMiller fc Scald 1 19791 ). the fourth the BH to total mass ratio, the 



fifth the BH mass in solar units, the sixth the primordial binary fraction. The seventh entry, (Am) Gog , is the 
snapshot-time-averaged value for the difference between mean main- sequen ce mass at the c enter of the cluster 
and mean main-sequence mass around r/j m , according to the definition of iGill et al.l (120081 ) . The eight entry, 
(Am) P09 , is the same quantity but normalized as discussed in section 2.2. Its standard deviation, minimum and 
maximum values are given in the last three colum ns. Snapshots fro m t = 7t r h to 9 t r h are used in the analysis. 
The first 12 entries are the simulations discussed in lGill et al.l (120081 ). The last 4 entries are additional runs, that 
are "randomized" clones of 16kmbh (generated with a different random number seed). 



